TT- and p-mesons, and their diquark partners, from a contact interaction 
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We present a unified Dyson-Schwinger equation treatment of static and electromagnetic properties 
of pseudoscalar and vector mesons, and scalar and axial-vector diquark correlations, based upon a 
vector-vector contact-interaction. A basic motivation for this study is the need to document a 
comparison between the electromagnetic form factors of mesons and those diquarks which play 
a material role in nucleon structure. This is an important step toward a unified description of 
meson and baryon form factors based on a single interaction. A notable result, therefore, is the 
large degree of similarity between related meson and diquark form factors. The simplicity of the 
interaction enables computation of the form factors at arbitrarily large spacelike Q^, which enables 



us to expose a zero in the p- meson electric form factor at Zq « ^GrUp. Notably, rpZ^ 
where rp, ro are, respectively, the electric radii of the p-meson and deuteron. 
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I. INTRODUCTION 

In numerous respects, tt- and p mesons are the simplest 
bound-states to study in QCD. That is, of course, suppos- 
ing that the framework employed is Poincare-covariant, 
capable of simultaneously describing light-quark confine- 
ment and dynamical chiral symmetry breaking (DCSB), 
and admits a symmetry-preserving truncation scheme. 
All these features are required because, amongst many 
other things, the pion is the lightest hadron and QCD's 
Goldstone mode, the p-meson couples strongly to two 
pions and is an important part of the photon's vac- 
uum polarisation, and modern facilities probe hadrons 
with momentum transfers far in excess of any reasonable 
constituent-quark- like mass-scale. 

The Dyson-Schwinger equations (DSEs) P, 01 provide 
an approach to hadron physics that is distinguished by 
its ability to satisfy these demands; and there is a large 
body of research that addresses tt- and p-meson prop- 
erties. For example, the analysis of static properties is 
reported in Refs. [3l-[23| and of interactions in Refs. [23 - 
|43[. There is nevertheless a need to return to this theme; 
namely, a programme aimed at charting the interaction 
between light-quarks by explicating the impact of dif- 
fering assumptions about the behaviour of the Bethe- 
Salpeter kernel on hadron elastic and transition form fac- 
tors nil. 

To expose the connection we remark that in quantum 
field theory a baryon appears as a pole in a six-point 
quark Green function. The pole's residue is proportional 
to the baryon's Faddeev amplitude, which is obtained 
from a Poincare covariant Faddeev equation that sums 
all possible quantum field theoretical exchanges and in- 
teractions that can take place between three dressed- 



quarks. A tractable truncation of the Faddeev equa- 
tion is based [i^ on the observation that an interaction 
which describes mesons also generates diquark correla- 
tions in the colour- 3 channel [5|. The dominant corre- 
lations for ground state octet and decuplet baryons are 
scalar (0'*') and axial-vector (1+) diquarks because, for 
example, the associated mass-scales are smaller than the 
baryons' masses and their parity matches that of these 
baryons. This is elucidated in Ref. [46j . 

At leading-order in a symmetry preserving truncation 
of the DSEs [111, [l^ , simple changes in the equations de- 
scribing TT- and p mesons yield expressions that provide 
detailed information about the scalar and axial- vector di- 
quarks; e.g., their masses H,S [13, El, 13113' ^^'^ electro- 
magnetic elastic (STj and transition form factors, which 
arc critical elements in the computation of a baryon's kin- 
dred properties. It is therefore natural to elucidate con- 
currently the properties of tt- and p-mesons and those of 
the scalar and axial- vector diquark correlations because 
it opens the way to a unified, symmetry-preserving ex- 
planation of meson and baryon properties as they are 
predicted by a single interaction. The potential of this 
approach is apparent in Refs. [3 El] but it has yet to 
be fully realised. For the present the best connection is 
provided by the less rigorous approach of Ref. [sl] , which 
uses more parameters to express features of QCD but 
also predicts and describes simultaneously a larger array 
of phenomena (5ll - [53| . 

Herein, as part of the programme outlined above, 
we describe results for a range of static and dynamic 
properties of these simplest u/d-mesons and -diquark- 
correlations as produced by a vector-vector current- 
current interaction that is mediated by a momentum- 
independent boson propagator; i.e., by the symmetry- 
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preserving regularisation of a contact interaction. Given 
the large body of work based on QCD-like vector-boson 
propagation that is already available, this study will pro- 
vide numerous points for comparison and contrast that 
are relevant to existing and planned experiments. 

In Sec.|TT] we describe a symmetry-preserving regular- 
isation and DSE-formulation of the contact interaction, 
following Refs. [l^, E^, EB] ■ Our scheme is such that con- 
finement is manifest, and chiral symmetry and the pat- 
tern by which it is broken are veraciously represented. 
In addition to the current-quark mass, the model has 
two parameters. In Scc. lIIII we describe results for tt- 
and p-meson electromagnetic elastic and transition form 
factors, computed using the rainbow-ladder truncation 
of the DSEs; with the analogous discussion of diquark 
correlations reported in Scc lIVl Section |V] provides a 
summary and perspective. 



II. CONTACT VECTOR-CURRENT-CURRENT 
INTERACTION 

A. Gap equation 

The typical starting point for a DSE study of hadron 
phenomena is the dressed-quark propagator, which is ob- 
tained from the gap equation: 



S{p) ^ — ij ■ p + m 



g Di,^{p- q)—-/^S{q)—T^{q,p), (1) 



wherein m is the Lagrangian current-quark mass, D^i, is 
the vector-boson propagator and is the quark-vector- 
boson vertex. Much is now known about D^^ in QCD 
[s^ - IHrj and nonperturbative information is accumulating 

on UllliMi- 

However, our goal is to build a stock of material that 
can be used to identify unambiguous signals in experi- 
ment for the pointwise behaviour of: the interaction be- 
tween light-quarks; the light-quarks' mass-function; and 
other similar quantities. Whilst these are particular qual- 
ities, taken together they will enable a characterisation 
of the nonperturbative behaviour of the theory underly- 
ing strong interaction phenomena [4l|, ■ We therefore 
elucidate predictions following from the assumptiorQ 

g'^D^y{p- q) = 6^y-^ , (2) 

where mc is a gluon mass-scale, and proceed by embed- 
ding this interaction in a rainbow-ladder truncation of 



^ This choice is the antithetical complement to that proposed in 
Ref. l62l : i.e., a i5- function in four- momentum space, which is con- 
fining because it provides a strong interaction that is independent 
of separation, . 



the DSEs, which is the leading-order in the most widely 
used, symmetry-preserving truncation scheme [isj . This 
means 



(3) 



in Eq. ([T]) and in the subsequent construction of the 
Bethe-Salpeter kernels. 

With this kernel the gap equation becomes 
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5-i(p) = z7-p + m+--^y^7.^(<z)7., 



(4) 



an equation in which the integral possesses a quadratic 
divergence, even in the chiral limit. If the divergence 
is regularised in a Poincare covariant manner, then the 
solution is 

S{p)-' =i-l-p + M , (5) 
where M is momentum-independent and determined by 
M f°° . 1 



M 



ds s ■ 



Af2 



(6) 



One must specify a reg ularisation procedure in order 
to proceed. We write [6j] 



M2 



JTIP ' 

where rir.uv are, respectively, infrared and ultraviolet reg- 
ulators. It is apparent from Eq. ([5]) that a nonzero value 
of =: l/Air implements confinement by ensuring the 
absence of quark production thresholds [gj, [11] . We note 
that since Eq. ^ does not define a renormalisable the- 
ory, Auv := I/tuv cannot be removed but instead plays 
a dynamical role and sets the scale of all dimensioned 
quantities. The gap equation can now be written 



M 



where C'"(A/2)/M2 



M 



(9) 



with r(a, y) being the incomplete gamma-function. 



B. Point-meson Bethe-Salpeter equation 

In rainbow-ladder truncation, with the interaction in 
Eq. © , the homogeneous Bethe-Salpeter equation for a 
colour-singlet meson is 



r(fc;P) 
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3 771q 



(2^ 



7mX(9;^)7/., (10) 
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where x{q; P) = S{q + P)r{q; P)S{q) and T{q; P) is the 
meson's Bethe-Salpeter amphtude. Smce the integrand 
does not depend on the external relative-momentum, 
k, then a symmetry-preserving regularisation of Eq. (jlOp 
will yield solutions that are independent of fc. It follows 
that if the interaction in Eq. ([2]) produces bound states, 
then the relative momentum between the bound-state's 
constituents can assume any value with equal probability. 
This is the defining characteristic of a pointlike composite 
particle. 

With a dependence on the relative momentum for- 
bidden by the interaction, the pseudoscalar and vector 
Bethe-Salpeter amplitudes take the general forrr0 [6^ 

r^(P) =i75i?.(/') + ^757-m(P), (11) 
^f.iP) - jJ^EpiP) + l^a^^P^FpiP) , (12) 
where P^jJ^ = and 7^ + 7^ = 7^- We observe that 



(13) 



However, it should be borne in mind that this is an arte- 
fact of the rainbow-ladder truncation; viz., even using 
Eq. ([2]), Fp[P) 7^ in any symmetry-preserving trunca- 
tion that goes beyond this leading-order [31 . We will 
see that the accident expressed in Eq. (|13p has material 
consequences. 



C. Ward-Takahashi identities 

No study of TT- or p-meson observables is meaning- 
ful unless it ensures expressly that the vector and axial- 
vector Ward-Takahashi identities are satisfied. The m = 
axial-vector identity states 



p^r5^(fc+, fc) = s-\k+)ij5 + i-f5S-\k) , 



(14) 



where Tr,^[k+, fc) is the axial- vector vertex, which is de- 
termined by 
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7^-T4 70X5^(9+, 9)70 



(15) 

We must therefore implement a regularisation that 
maintains Eq. (|14p . This requirement is readily found 
to entail the following two chiral limit identities (40j : 



8 M 
6 mi 



d^q 



1 



1 



AP 



= 



(2^ 



P-q+ 



P-q 



q\ + AI^ q'^ + AP 



(16) 
(17) 



^ We assume isospin symmetry throughout and hence do not in- 
clude the Pauli isospin matrices exphcitly. 



which must be satisfied after regularisation. Analysing 
the integrands using a Feynman parametrisation, one ar- 
rives at the follow identities for = = to: 



, r 16 M 

6 TOq 



d'^q 



1 







(27r)4 [q2 + AP] ' 

dS ^q^ + AP 
(27r)4 [g2 + A/2]2- 



(18) 
(19) 



Equation (|18p is just the chiral-limit gap equation. 
Hence it requires nothing new of the regularisation 
scheme. On the other hand, Eq. (fT9|) states that the 
axial-vector Ward-Takahashi identity is satisfied if, and 
only if, the model is regularised so as to ensure there 
are no quadratic or logarithmic divergences. Unsurpris- 
ingly, these are the just the circumstances under which a 
shift in integration variables is permitted, an operation 
required in order to prove Eq. (jl4p . 

It is notable, too, that Eq. is valid for arbi- 
trary P. In fact its corollary, Eq. ([T5| . may be used 
to demonstrate that in the chiral limit the two-fiavour 
scalar-meson rainbow-ladder truncation of the contact- 
interaction DS Es p roduces a bound-state with mass 
TOct ~ 2 Ad (46l |67{ . In the presence of a momentum- 
dependent dressed-quark mass function, one could re- 
verse this association and define a chiral-limit dressed- 
quark constituent-mass as one-half the mass of the light- 
est rainbow-ladder scalar meson. This procedure yields 
M° ~ 0.3 GeV, as may readily be determined from 
Ref. 0. 

The second corollary, Eq. P7)) . entails 

= / da [C'"(a;(M2,a,F2)) + C\^{uj{AP ,a, P^))] , 

(20) 



with 



(21) 



uj{AP,a, p2) = AP + a{l - a)P^ , 
Ci"(z) = -z(d/dz)C*"(z) 

= z [r(0,MV2j - r{0,APrl)] . (22) 

The vector Ward-Takahashi identity 

P^iri{k+, fc) = S-\k+) - S-\k) , (23) 

wherein VJ^ is the dressed-quark-photon vertex, is crucial 
for a sensible study of electromagnetic form factors [l^ . 
Ideally, the vertex needs to be dressed at a level con- 
sistent with the truncation used to compute the bound- 
state's Bethe-Salpeter amplitude [2§|. In our case this 
means the vertex should be determined from the follow- 
ing inhomogeneous Bethe-Salpeter equation: 



r,,(Q)-7^. 
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3 TOq 



dS 

(27r)4 



laXf^iq+.qha, (24) 



where x^((?+,q) = S{q + P)r^(Q)S'(g). Owing to the 
momentum-independent nature of the interaction kernel, 
the general form of the solution is 



r^{Q) = ij:PT{Q') + if:PL{Q').. 



(25) 
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FIG. 1. (Color online) Dressing function for the transverse 
piece of the quark-photon vertex; viz., Pt{Q^) in Eq. (|27|l . 



where Q^i^^ = and 7 J + 7^ = 7^- This simplicity 
doesn't survive with a more sophisticated interaction nor 
with Eq. ^ beyond rainbow-ladder truncation [i^ . 
Inserting Eq. ([25]) into Eq. , one readily obtains 



Pl{Q^)^1, (26) 
owing to Eq. (|17p . Using this same identity, one finds 

1 



with (Ci(z) =Ci(z)/z) 



l + i^^(Q2)' 



(27) 



37r2m^ 



X/ daa(l-a)Q2Ci (w(M2,a,g^)). (28) 
/o 



Plainly, 



Pt(Q' = 0) = 1. 



(29) 



so that at — in the rainbow-ladder treatment of the 
interaction in Eq. ^ the dresscd-quark-photon vertex is 
equal to the bare vcrtexH 

However, this is not true for ^ 0. In fact the trans- 
verse part of the dressed-quark-photon vertex will display 
a pole at that < q for which 



(30) 



This is just the model's Bethe-Salpeter equation for the 
ground-state vector meson. 

In Fig. [T] we depict the function that dresses the trans- 
verse part of the quark-photon vertex. The pole as- 
sociated with the ground-state vector meson is clear. 



^ Equations II26I I. I I29I I guarantee a massless photon and show 
that our regularisation also ensures preservation of the Ward- 
Takahashi identity for the photon vacuum polarisation |68|| . 



Another important feature is the behaviour at large 
spacelike- Q 2; namely, Pt{Q^) — >■ 1~ as oo. This 

is the statement that a dressed-quark is pointlike to a 
large-Q^ probe. The same is true in QCD, up to the loga- 
rithmic corrections which are characteristic of an asymp- 
totically free theory (29l |. 



D. Bethe-Salpeter kernels for tt and p 

At this point we can write the explicit form of Eq. (flO 
for the pion: 
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JCef 










K-FE 


K-FF 







, (31) 



where 

K^EE 



da [C''iu}{M^,a,~ml)) 



JCef = ^fn 
1 



+2a{l-a)mlC,{ij{M^a,~ml)) , (32) 
1 

2 „ _2 ^ 



daCi {u!{M , a, —m^)), 



(33) 



JCfe 

K-FF 



^M^ daC'^{Lo{M\a,~ml)), (34) 



= -2/C 



EE ■ 



(35) 



This is an eigenvalue problem for the pion mass-squared, 
ml. NB. We used Eq. ^ to arrive at Eq. (|35|) . 

The explicit form of Eq. pU)) for the p-meson, whose 
solution yields its mass-squared, is 

l + /v^(-m2)=0, (36) 

where is given in Eq. (pS)) . 

In the computation of obscrvablcs, one must use the 
canonically-normalised Bethe-Salpeter amplitudes. For 
the rainbow-ladder pion this means that Ftt is rescaled 
to satisfy 

Pf^ = NctTj ^T.(-P)^S{q + P) r.(P) S{q) , 

(37) 

which, in the chiral limit, becomes 



1 



Nr 1 



4^2 C,{M^-tI, tI,)E,[E, - 2F^]. (38) 

For the rainbow-ladder p-meson, on the other hand, the 
vector meson analogue of Eq. (|37|) requires that 



1 



^ dz ^ 



(39) 



In terms of the canonically normalised Bethe-Salpeter 
amplitudes, the leptonic decay constants of the tt- and 
p-mesons are respectively given by 

1 3 p2 



fp 



^ Pp T.^ I 2\ 



(40) 
(41) 
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TABLE I. Results obtained with (in GeV) mc = 0.132 , Air = 
0.24 , Auv — 0.905, which yield a root-mean-square relative- 
error of 13% in comparison with our specified goals for the 
observables. Dimensioned quantities are listed in GeV. 



m 


-Ejr -Fit Ep 


M kJ'^ mp fp 




0.007 


3.568 0.459 1.520 
3.639 0.481 1.531 


0.358 0.241 0.919 0.100 0.130 
0.368 0.243 0.140 0.928 0.101 0.129 



Another important low-energy property is the in-pion 
condensate^ 



F^JC 



EE 



(42) 



In the chiral limit k-,^ — > — ~{qq); i.e., the so-called 
vacuum quark condensate j69| . Moreover, in this limit, 
too, one can readily verify that [i^ 



m=0 M 



(43) 



which is a particular case of one of the Goldberger- 
Treiman relations proved in Ref. [l^, and Ftj{P = 0) 
satisfies a similar identity. 



III. TT AND p ELASTIC AND TRANSITION 
FORM FACTORS 



In order to compute the form factors we need to fix the 
model's two parameters; namely, ma and A^vE We do 
that by performing a least-squares fit in the chiral limit 
to MO = 0.40 GeV, k° = (0.22 GcV)^, /o = 0.088 GeV, 
m° = 0.78 GeV and /" ^ 0.15 GeV. This procedure 



yields the results in Table HI 



A. TT-meson elastic form factors 

We are solving the interaction of Eq. ([2]) in the 
rainbow- ladder truncation; i.e., at leading-order in 
the nonperturbative symmetry-preserving truncation of 
Refs. [U, [l^l. At this order the generalised impulse ap- 
proximation is computed for three-point scattering pro- 
cesses [IHl, such as elastic form factors. An analy- 
sis of the associated triangle diagram yields the formu- 
lae in Sec. lA II and the computed result is depicted in 
Fig. [2] Two features are immediately apparent; viz., the 
pole associated with the p-meson at timelike momentum. 



E 




FIG. 2. (Color online) F°™((3^) computed in rainbow-ladder 
truncation from the interaction in Eq. solid curve - fully 
consistent, i.e., with a dressed-quark-photon vertex so that 
the p-pole appears; and dashed curve - computed using a 
bare quark-photon vertex. Dotted curve ~ fit to the result 
in Ref. [sJl , which also included a consistently-dressed quark- 
photon vertex and serves to illustrate the trend of contempo- 
rary data. 



which is a consequence of dressing the quark-photon ver- 
tex; and a momentum-independent interaction produces 
F-k{Q^) = constant as — )■ oo. The following function 
is a valid interpolation of the full result on the domain 
shown: 

pcm.^2^ interpolation 1 + 0.33 + 0.024 

" ' 1 + 1.20 Q2 + 0.053 Q4 ^^^> 

In Tablellll we report the pion charge radius: 
d 



= -6 



F^iQ') 



(45) 



The result is less than experiment (r^r = 0.672 ± 0.008 fm 
[toII). This owes in small part to our omission of pseu- 
doscalar meson rescattering effects [7l| but more to the 
large value we obtain for the p-meson's mass. It cannot 
be remedied in our symmetry-preserving rainbow-ladder 
treatment of Eq. ^ because all dimensioned quantities 
are too closely tied to the value of M. An interaction 
which preserves the one-loop renormalisation group be- 
haviour of QCD provides decoupling between the 
values of ultraviolet and infrared phenomena, such as 
and K^. 



B. /o- meson elastic form factors 



* There is an analogous in-p-meson condensate but that will be 

discussed elsewhere. 
5 We fix Air = 0.24 GeV Aqcd since rqcD := 1/Aqcd ~ 0.8 fm 

is a length-scale typical of confinement. 



The J '-^ = 1 p-meson has three elastic form factors 
and we follow Ref. [1^ in defining them. Denoting the 
incoming photon momentum by Q, and the incoming and 
outgoing p- meson momenta by = K — Q/2 and p-^ = 
K + Q/2, then K ■ Q = 0, + = -m^ and the 
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TABLE II. Row 1: Form factor radii (in fm), and magnetic 
and quadrupole moments for the p-meson, GIAQ'^ = 0) and 
Gq{Q^ — 0) respectively, computed with (in GeV) m = 0.007, 
niG = 0.132, Air = 0.24, A^v = 0.905. For a structure- 
less vector meson, fj. = 2 and Q = —1 [T^I. The next 
four rows list results reported elsewhere. Experimentally, 
= 0.672 ±0.008 fm [t^]. (NB. None of the quoted computa- 
tions included contributions from nonresonant pseudoscalar- 
meson final-state interactions and hence agreement with the 
experimental value of should be seen as a defect of the as- 
sociated model [tJ]. The nature of this flaw is understood 
within the DSE context, where such contributions can be 
viewed as computable corrections to the rainbow-ladder trun- 
cation d^.) The last two lines report results for the scalar 
and axial-vector diquark correlations. Here the magnetic and 
quadrupole moments should be multiplied by the relevant 



charge factor; viz., ej^u) 



4 
3 ' 



'-{ud} 



— i and e 



{dd} 



This work 


0.45 


' p 

0.56 


' p 

0.51 


' p 

0.51 


Mp 
2.11 


Qp 
-0.85 


Ref. [39] 


0.66 


0.73 






2.01 


-0.41 


Refs. [73, 74] 


0.56 


0.61 






2.69 


-0.84 


Refs. [75, 76] 


0.66 


0.61 






2.14 


-0.79 


Refs. [77, 78} 


0.66 


0.52 






1.92 


-0.43 


This work 


'■0+ 
0.49 


^1 + 

0.55 


0.51 


0.51 


Mi+ 
2.13 


Qi+ 

-0.81 


Ref. -[37] 


0.71 













p-^-p vertex can be expressed: 

3 



^2m2 
Pi 



2ml 



2ml 



(46) 
(47) 

(48) 

(49) 

where 'P'^^{p) = S^i^ — p^p^ /p^ . A symmetry-preserving 
regularisation scheme is essential here so that the follow- 
ing Ward-Takahashi identities are preserved throughout 
the analysis: 

Qx^xMK,Q)=Q (50) 

p^Aa,^.(a-, g) = = pIa^mk. Q) ■ (51) 

The electric, magnetic and quadrupole form factors are 
constructed as follows: 



(52) 

GM(Q') = -i^2(Q'), (53) 
Gq(Q2) = Fi(Q2) + F2{Q^) + [1 + 77] F^{Q^) , (54) 



where r] = Q'^/[Am1\. In the limit _j. g, these form 
factors define the charge, and magnetic and quadrupole 
moments of the p- meson; viz., 



G^MiQ' 



0) = 1 , (55) 
0)=/.p, G^(Q2 = 0) = gp. (56) 

It is readily seen that Eq. (|55p is a symmetry constraint. 
One has GsiQ^ =0) = Fi{Q^ = 0) and 



A{K, Q)''^'' 2K, Vl^ (K) (K) F, (0) . 
Using Eqs. (|JTO| . this becomes 



(57) 



(0y*7. ^^Sie + K)z^^ S(i) . (58) 



The right-hand-side (rhs) is simply the analogue of 
Eq. p7p for the rainbow-ladder vector meson. Hence, 
when Ep is normalised according to Eq. ([39|) and so 
long as one employs a symmetry-preserving regularisa- 
tion procedure, the rhs is equal to K\V^^{K) and thus 
-Fi(0) = l. 

We compute the form factors using the formulae in 
Sec. lA 21 In Table |lT] we report form factor radii, and 
the magnetic and quadrupole moments. The comments 
following Eq. (|45|) are also relevant to the magnitudes of 
the yO- meson radii. An interpretation of the ratio r^/r^ = 
0.80 determined from the Table is complicated by the 
fact that we have consistently used the rainbow-ladder 
truncation; but in this case alone Fp{P) = 0, whereas 
Ftt{P) 7^ always and Fp{P) 7^ in all other truncations. 
We observe therefore that r^r = 0.51 fm if one artificially 
sets Ft^{P) = 0, in which case r^Jr-p = 0.92. More- 
over, the DSE computation in Ref. [391 1 . which employs a 
QCD-based interaction, produces r^^/rp = 0.90; and in 
combination, the more phenomenological DSE studies of 
Refs. [7i,[l3l yield r-^/vp = 0.92. 

Our computed p-meson electric form factor is plotted 
in Fig. [31 It displays a zero at = 5.0 GeV^ and re- 
mains negative thereafter. Given that the deuteron is 
a weakly-bound J = 1 system, constituted from two 
fermions, and its electric form factor possesses a zero [79| . 
it is unsurprising that C^i^Q^) exhibits a zero. It is no- 
table in addition that the deuteron's zero is located at 
zS := VQ^ = 0.8 GeV, so that 



p E 



(59) 



where td is the deuteron's radius. An interpolation valid 

on Q2 g [-777,2, 10 GcV^] jg 



G%{Q^) 



2\ interpolation 



1 - 0.20 q2 



1 + 1.15 Q2 _ 0.013 Q4 



(60) 



In Fig.[3lwe also depict the magnetic and quadrupole 
form factors of the p-meson, both normalised by their 
values = 0. Notably, neither of these two form fac- 



tors change sign: for Q > 



G^iiQ ) is positive 
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Q, 




2 4 6 8 
Q2 (GeV^) 
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FIG. 3. (Color online) Solid curve - p-meson electric form 
factor, G'^{Q^), which exhibits a zero at Q2 = 5.0 GeV^ 
(It is notable that 1 — |?7 = for = Qm^ = 

5.2 GeV^) The dashed curve, G'Ij{Q^) / ^.p, and dot-dashed 
curve, Gq{Q'^) / Qp, are almost indistinguishable, as empha- 
sised by the dotted curve, [G5f(Q^)/Mp]/[G{^(Q^)/Qpl- The 
charge radii, and magnetic and quadrupole moments are given 
in Table [m NB. All form factors exhibit a pole at = — rUp 
because the quark-photon vertex is dressed as described in 
Sec.HTCl 



definite and G'^iQ'^) is negative definite. Furthermore, 
over this entire domain of Q^, these form factors exhibit 
a very similar Q^-dependence, which is made especially 
apparent via the dotted-curve in Fig. [31 Interpolations 
vahd on e [-ml, lOGeV^] are 



i^P interpolation 



2.11 + 0.021 Q2 
H^1.15Q2_ 0.015 Q4 
0.85 + 0.038 Q2 



ip //^2\ interpolation 

''^^^ ' ^ ^l + 1.17Q2 + 0.014Q4 



(61) 
(62) 



The similar momentum-dependence of G^^ and Gq re- 
calls a prediction in Ref. [t^ ; namely, 

Ge{Q') ■■ Gm{Q') : Gq(Q2) 1 - : 2 : -1 (63) 

in theories with a vector- vector interaction mediated via 
bosons propagating as l/fc^ at large-fc^. Our computed 
ratio rjv//Q := G''j^j{Q^) /Gq{Q'^) conforms approximately 
with this prediction on a large domain of Q^; e.g., 





rM/Q -2-4 



10 10^ 10^ 
-2.54 -2.38 -2.17 



(64) 



However, at = lO-^GeV^, tm/q = -1.28. Moreover, 
the remaining two ratios are always in conflict with the 
prediction; and closer inspection reveals that even the 
apparent agreement for G'^j{Q'^)/Gq{Q'^) is accidental, 
since Eqs. (j63p are true if, and only if, 

Fi{Q') : F2(Q2) : Q^F^iQ^) '^'=°° 1 : -2 : ; (65) 

and none of these predictions are satisfied in our compu- 
tation. 





1.0 




0.8 




0.6 




0.4 




0.2 




0.0 




2 4 6 8 
(f (GeV^) 



10 



FIG. 4. (Color online) Solid curve - the full result for 
G''^''(Q^); and dashed curve - G"'^''(0^) obtained with 
_FV(P) = 0. Experimentally [t^I, the partial width for 
— j> 7r"'"7 is 68 ± 7keV, which corresponds to [sj] g^-yp ~ 
(0.74 + 0.05) mp. This is in fair agreement with our computed 



The mismatch originates, of course, with Eq. ^ and 
the concomitant need for a regularisation procedure in 
which the ultraviolet cutoff plays a dynamical role. If one 
carefully removes Auv — ^ oo, Eqs. (|65p are recovered but 
at the cost of a logarithmic divergence in the individual 
form factors. We conclude therefore that a vector-vector 
contact interaction cannot reasonably be regularised in a 
manner consistent with Eq. (|63p . 

In closing this subsection we reiterate that it is only in 
the rainbow-ladder truncation that Fp(P) = 0. Therefore 
in connection with the p-meson's form factors, material 
changes should be anticipated when proceeding beyond 
this leading-order truncation. 



C. p-TT transition form factor 

This transition is closely related to the 7*7r7 transition 
form factor, whose behaviour in connection with Eq. ([2]) 
was analysed in Ref. [1^. The interaction vertex is ex- 
pressed in Eq. (|A19[) and defines a single form factor; viz., 



r;j''(fci,fc2) 



.fsha.k2f3G^-"'iQ^), (66) 



where kf — Q , k2 — ^tTip. The coupling constant, gn-^p: 
is defined such that G'^^p{Q^ — 0) = 1; and explicit 
formulae for computing this form factor are provided in 
App.El 

Our computed form factor is depicted in Fig.|4l Nat- 
urally, because the quark-photon vertex is dressed (see 
Fig.[D), the transition form factor exhibits a pole at 
~ —m^p, which we have not displayed. An interpola- 
tion vahd on e [-nri^p, 10 GeV^] is 

r-<T^lP(rfi\ interpolation 1 + 0.37 + 0.024Q^ 

^ " 1 + 1.29 Q2 + 0.015 Q4- ^^'^ 
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FIG. 5. (Color online) Solid curve - full result for scalar- 
diquark elastic electromagnetic form factor; and dashed curve 
- result obtained without dressing the quark-photon vertex. 
The computed mass of the diquark is rriqq^^ = 0.776 GeV and 
the charge radius is given in Table [III 



In the neighbourhood of = 0, the form factor is 
characterised by a radius-like length-scale; viz., 



TTfp 



-6 



(0.46 fm)^ 



(68) 



which is almost indistinguishable from both r^r = 0.45 fm 
in Table [n] and the anomaly interaction radius defined in 
Ref. m 

'. viz., T*Q — 0.48 fm. On the other hand 



lim G"''^''(Q2) ^ 0.11 



(69) 



owing to the presence of the pion's pseudovector com- 
ponent, a result in keeping with the pointlike nature of 
bound-states generated by a contact-interaction [i^ \^ . 



IV. 0+- AND 1+-DIQUARK ELASTIC AND 
TRANSITION FORM FACTORS 

A. Scalar-diquark elastic form factor 



CD 



CD 

6 

CD 




FIG. 6. (Color online) Solid curve - Pseudovector-diquark 
electric form factor, {Q^), which exhibits a zero at = 
6.5 GeV^ (In this case 1- §77 = for = 6m^+ = 6.7GeV^ 
given the computed mass of 1.06 GeV.) The dashed curve, 

G\i (Q^) / , and dot-dashed curve, Gq^(Q^)/Qi+, are al- 
most indistinguishable, as emphasised by the dotted curve, 

[Gm iQ^)/l^p]/[GQ iQ^)/Qi+]- The charge radu, and mag- 
netic and quadrupole moments are given in Table [III NB. 
All form factors exhibit a pole at — — rrip because the 
quark-photon vertex is dressed as described in Sec. Ill Cl 



Our result for the scalar diquark elastic electromag- 
netic form factor is presented in Fig.[S] An interpolation 



valid on e [-ml, lOGeV^] is 



TTiom 
^0+ 



, „2\ interpolation 1 1 -I- 0.25 



0.027 Q4 



3 1-f 1.27Q2 + 0.13Q4 



(72) 



The normalisation is different but the momentum- 
dependence is similar to that of F°'°^. This is indicated, 
too, by the ratio of charge radii; viz., rQ+frT^ = 1.08, 
which may be compared to the value of 1.09 obtained 
in Ref. [33| and contrasted with the value of 0.8 in [soj . 
In the absence of the scalar-diquark Bethe-Salpeter am- 



plitude's pseudovector component, 



0.51 fm; i.e., an increase of 6%. 



1%+ 



0, we find 



In the context of the interaction in Eq. ([2]), a detailed 
discussion of the relationship between pseudoscalar- and 
vector-mesons and scalar- and axial-vcctor-diquark cor- 
relations may be found in Ref. [46| . Using the informa- 
tion provided therein, it is straightforward to show that 
in rainbow-ladder truncation the electromagnetic form 
factor of a scalar diquark is readily obtained from the 
expression for F°™(Q^). Namely, 



""0+ 



, (70) 



where the scalar-diquark Bcthe-Salpeter amplitude is ex- 
pressed via (C = 7274 is the charge-conjugation matrix) 



r,,„+ iP)C^ - 75 {P) + ^7 • PF,,^,^ (P) 



(71) 



B. Pseudovector-diquark elastic form factors 

From the above observations it will be apparent that 
the rainbow-ladder results for the {ud} axial- vector di- 
quark elastic form factors may be obtained directly from 
those of the p-meson through the substitutions 



1 



(73) 



The momentum-dependence of the form factors for the 
{uu} and {dd} correlations is identical but in these cases 
the normalisations are, respectively, | and — |. 

We depict the axial-vector diquark form factors in 
Fig. El They are similar to but distinguishable from those 
of the p-meson, falling-off a little less rapidly owing to the 
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FIG. 7. (Color online) Solid curve - momentum-dependence 
of full result for axial-vector-scalar-diquark transition form 

factor, G°^^^^(Q^); and dashed curve - result for G''^''(Q^) 
in Fig.m The different rates of evolution are typical of me- 
son cf. diquark form factors computed herein. Note that 
e{ud}3o+7i+"T-g9i+ = e{„tj}0.74 = 0.25. 



larger mass of the axial-vector diquark. Interpolations 
vahd on e [-"ip, lOGeV^] are 

/ interpolation 1 0.16 Q'^ 

{Q ) = 



^2^+ ,^2\ interpolation 
[Q ) = 

^]^+ ,^2\ interpolation 



1 + 1.1702 + 0.012 Q4 
2.13-0.19Q2 

1 + 1.07 02 _ 0.10 g4 ' 

0.81 - 0.029 g2 



(74) 
(75) 
(76) 



1-M.lig2_ 0.054 Q4 

from which the particular pseudovector diquark form fac- 
tors arc obtained after multiplication by the appropriate 
charge factors, listed in Table HIl 



C. l~'~-0~'~ diquark transition form factor 

Owing to the flavour structure of the scalar diquark, 
this transition can only involve the {ud} axial- vector di- 
quark. It is described by a single form factor, which can 
be introduced through 



1 go+7i+ 

3 TO, 



(77) 

and one may readily determine that in rainbow-ladder 
truncation 

G°^''i^(Q2) 

Computation of the form factor is straightforward and 
the result is depicted in Fig.[71 An interpolation valid on 

g2 e [-TO2,10GeV^] is 

^0+71+ /^2\ interpellation 1 + 0.10 , , 

> - 1 + 1.073 02- ^'^> 



The associated transition radius is 

ro+^i+ = 0.48fm, (80) 
which is 5% larger than r-^-^p in Eq. ()68p , and 

lim G°^^i^(Q2) = 0.049, (81) 
just under one-half of the value in Eq. (|69p . 

V. EPILOGUE 

We described a unified Dyson-Schwinger equation 
treatment of static and electromagnetic properties of 
pseudoscalar and vector mesons, and scalar and axial- 
vector diquark correlations based upon a vector-vector 
contact-interaction. Isospin symmetry was assumed, 
with TO„ = TOrf = TO = 7McV producing a physical pion 
mass; and two parameters were used to define the gap- 
and Bethe-Salpeter. In a comparison with relevant static 
quantities, we recorded a value of 13% for the overall 
root-mean-square relative-error. 

A basic motivation for our study is the need to doc- 
ument a comparison between the electromagnetic form 
factors of mesons and those diquarks which play a mate- 
rial role in nucleon structure because this is an important 
step toward a unified description of meson and baryon 
form factors based on a single interaction. A notable 
feature of our results, therefore, is the large degree of 
similarity between related form factors. For example, we 
find that it would be a good practical approximation to 
assume equality of related radii: ro+ ~ r.^ and ri+ w rp. 

As has previously been observed, a fully-consistent 
treatment of the contact interaction produces a pion elec- 
tromagnetic form factor that approaches a nonzero con- 
stant value at large spacelike momenta. On the other 
hand, owing to a peculiarity of the rainbow-ladder trun- 
cation, which prevents the appearance at this order of 
a tensor component for the p-mcson produced by a con- 
tact interaction, the p-meson form factors approach zero 
at large spacelike momenta. This accident means that a 
comparison with QCD-based DSE calculations can mean- 
ingfully be interpreted. In a comparison with the most 
sophisticated such study, the form factors produced by 
the contact interaction are harder although the ratio 
rjr/r^ is similar. Moreover, the contact interaction's 
simplicity allows one to readily compute the p-meson 
form factors at arbitrarily large spacelikc and ex- 



6m p. 



pose a zero in the electric form factor at 
Notably, tdZq ~ "i^p^qi where td and Zq arc, respec- 
tively, the deuteron's radius and the location of the 
zero in its electric form factor. The p-meson's magnetic 
and quadrupole form factors are positive- and negative- 
definite, respectively. We reiterate that the behaviour 
of all pseudovector-diquark form factors is semiquantita- 
tively the same. 

At the core of our analysis is a symmetry-preserving 
treatment of a vector-vector contact interaction. This 
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has now been used in the completely-consistent compu- 
tation of the hadron spectrum, and meson and diquark 
form factors. The foundation has thus been laid for the 
computation of baryon elastic and transition form fac- 
tors, which will provide information that is crucial for 
the use of experimental data on such observables as a 
tool for charting the nature of the quark-quark interac- 
tion at long-range (44 1. 



tion Program of the Chinese Academy of Sciences, Grant 
No. KJCX2.YW.W10. 



Appendix A: Form Factor Formulae 

This appendix is a repository for the formulae we have 
used to compute the form factors. 
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1. Elastic pion form factor 



+ e^f^t:"},aQ') + F'n-^FiQ')] ,(Ai) 

where Pt(Q^) is given in Eq. (j27p and 



rpem 



3 
4^ 

3 
2^ 



daC^ {Lu{M^,a,Q^)) + 2ml / dadfi aC^ {oJ2{M^ ,a, (i,Q^ ,ml)) 



daC^{uj{M',a,Q'))+ dadp a{aQ' - 2mi)C^ {uj2{M\a, P,Q\mi)) 



rpcm 



A{a, Q\ml)C, (u,2{M\a, /3, Q\ ml)) 
+ [B{M\ a, /3, Q\ ml) - A{a, Q\ ml) c.2(A/^ a, /3, Q\ml)] t^{u:2{.M\a, /3, Q^ )) 



(A2) 

(A3) 

(A4) 
(A5) 



with 



B{M\a,l3, Q\mi) = aQ^M' + M'mi{a - 2) + ami{aQ'[l -a~ 2/3(1 - /?) -f 3a/3(l - /?)] - (1 - (A6) 



^(a,g^m^) 



l„4(2_3a), (A7) 



2. Elastic p-meson form factors 



W2(M', a, /?, 0', to' ) = + g'a2/3(l - /?) 

-a{\~a)ml, (A8) 



In generalised impulse approximation the p-7 vertex in 
Eq. becomes 

--|-i?p( V)7J S{q + /) Pr(Q')m 



(27r4) 

^S{q + f)Ep(j>^)^lS{q)., 



(AlO) 



where C™{z) is defined after Eq. (O; 0^{z) and 
ij{M^,a, Q2) in Eqs.dUl), (1211); C'"(z) after Eq. 1^; and 



(A9) 



withC2"=C^"(z)/22. 



where Ep is the canonically-normalised p-mcson Bethe- 
Salpcter amplitude. Explicit expressions for the scalar 
functions Fi,2,3(Q') can be obtained via contraction with 
any three sensibly chosen projection operators; and the 
subsequent use of Feynman parametrisations yields 



(All) 



F,{Q^) = -^El j dadPa [ACT{i02) 

+ [Bi - AiUJ2]C'2{^^2) , 
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where = T,{M^,a,P,Q^,ml), = A^, Bf, viz., _ 2a(l - 2/3)(577t^ + Q^) 



Ai=2-a, (A12) -J • ■■-p- 

_ m2(a(10/? - 7) - 4) + Q^a[2p - 1) -a^l3{2 - a)(l - /3)Q2] , (A15) 

-42 ^ , (A13) 

L I 



TO^Bz = ml[M^{-A -7a + lOaP) - m^(-l + a)a{l - 7a - 6/3 + lOa/^)] 

+a[Af2(-l + 2/3) + mla{-l + 2/3 + a[l + l3- 5/3^ + 2I3^])Q^] , (A16) 
{4ml + = 4a [ml{5M^l ~ 2/3) + m2(-l + a)[3 - 6/3 + a(-5 - 6/3 + 16/?^)]) 

+ (Af2(l - 2/3) - m2a[-l + a - 2/3 + 3a/3 + 4/3^ - 7a;32 + 2al3^])Q'^] . (A17) 



3. Vector-pseudoscalar transition form factor 

The interaction vertex describing the n-p transition 

T;:jp{ki,k2) = ^ e^,,^pki^k2p G^"'P{Q^) (A18) 

TTlp 

= t^o/^^j T,{-P)S{e2) Pt{Q') lip 
xS{h2)iK{k2)S{ii), (A19) 

where the incoming p-meson has momentum k2 , the pho- 
ton has momentum ki = Q, the outgoing pion has mo- 
mentum P = (fci + fc2); and £i = £ — ki, £2 = £ + k2, 
£12 = £ — ki + k2. In this instance the kinematic con- 
straints are 



kl = -ml, kj^Q^, 2ki-k2 = ml-ml-Q^ . (A20) 



Given the structure of the pion's Bethe-Salpeter am- 
phtude, one may write 



G^^P{Q^) = Gl^'iQ^) + G'P'iQ^) , (A21) 



wherein 

GE'^iQ') = f dadpaC)^{oj,) , (A22) 

271"^ Jo 

07"^') = -^^ fdadPa [fr'CUu:,) 

J 

^{ir'-^^!r')c^{^^)\. (A23) 

with 

= -a [a/3(l - P)ml + (1 - a)(l - P)ml 
-(l-a)/3g2] (A24) 

and 

= 2 - 3a , (A25) 
/;^'' = (2-a)(M2+a2/3(l-/3)m2) 

+ {l-af{a{l-p)ml~al3Q^). (A26) 

The vertex in Eq. (jA19[) is intimately connected with 
the Abehan anomaly, which describes the process 7r° — > 
77 and associated transition form factors. The manner 
by which all aspects of the anomaly may be reproduced in 
the model we're considering is detailed in Sees. III. A and 
III.B.2 of Ref. [l^]. In the present context, consistency 
with the anomaly requires that in Egs. (jA22p . (jA23p . 
Sp = Ep/M-K-y-y, with M-ir-y-i defined such that G'"^'''{Q'^ = 
0) = 1/2, and Gy(Q'') ^ Gp^''{Q^) - Gp^''{0). Both 
modifications are necessary in order to correct for the dy- 
namical role played by the ultraviolet cutoff in a contact- 
interaction theory. 
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